Fractality of eroded coastlines of correlated landscapes 
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We investigate through numerical simulations of a simple sea-coast mechanical erosion model, the 
effect of spatial long-range correlations in the lithology of coastal landscapes on the fractal behavior 
of the corresponding coastlines. In the model, the resistance of a coast section to erosion depends 
on the local lithology configuration as well as on the number of neighboring sea sides. For weak 
sea forces, the sea is trapped by the coastline and the eroding process stops after some time. For 
strong sea forces erosion is perpetual. The transition between these two regimes takes place at a 
critical sea force, characterized by a fractal coastline front. For uncorrelated landscapes, we obtain, 
at the critical value, a fractal dimension D — 1.33, which is consistent with the dimension of the 
accessible external perimeter of the spanning cluster in two-dimensional percolation. For sea forces 
above the critical value, our results indicate that the coastline is self-affine and belongs to the KPZ 
universality class. In the case of landscapes generated with power-law spatial long-range correlations, 
the coastline fractal dimension changes continuously with the Hurst exponent H, decreasing from 
D — 1.34 to 1.04, for H = and 1, respectively. This non-universal behavior is compatible with the 
multitude of fractal dimensions found for real coastlines. 

PACS numbers: 89.75.Da,92.40.Gc,64.60.al,64.60.ah 



I. INTRODUCTION 

Since the introduction of the concept of fractals by 
Mandelbrot [l^ , scale invariant behavior has been identi- 
fied and investigated in many geological and geophysical 
phenomena [2], including the frequency distributions of 
earthquakes [3] and volcanic eruptions , the size distri- 
bution of rock fragments [51 |S] and mineral deposits [7] , 
and the topography of river networks [51 H] , rivers deltas 
[To] and rocky coastlines [11] [12] . More recently, the ge- 
ometry of watersheds, namely, the lines separating adja- 
cent drainage basins (catchments), has been reported to 
display also typical fractal features [TBHT5] . with impor- 
tant implications to water management [16', landslides 
[T71 IT8j , and flood prevention [19j . 

Erosion is certainly one of the most remarkable ex- 
amples of a geological process that naturally generates 
diverse self-similar structures. In particular, the erosion 
of a coastline by the sea [20] constitutes a rather rich 
phenomenon. Due to the action of underlying geological 
processes (e.g., tectonics and volcanic events), topologi- 
cal and lithological properties of coastal landscapes are 
generally heterogeneous as well as long-range correlated 
in space. As a consequence, the resistance to erosion 
must be considered as a spatially dependent parameter. 
One should therefore expect that the self-similar geome- 
try of coastlines should emerge from an intricate interplay 
between these landscape properties and the sea force. 



In a study by Sapoval et al. |21j . an erosion model 
is proposed to show how surviving coastlines can dy- 
namically evolve to self-similar objects by means of a 
self-organized critical process (22, ,23^ . By applying this 
model to spatially uncorrelated landscapes, they recov- 
ered, at the critical steady-state, a fractal dimension 
D ~ 1.33 that is frequently observed in real systems 
[T2] . This dimension is different from 7/4, the dimen- 
sion of the external perimeter of percolation |24j , but it 
is consistent with 4/3, the one for the accessible external 
perimeter |25j . This value, however, cannot be taken as 
universal, since, in fact, a multitude of fractal dimensions 
has been measured for real coastlines [26] . 

In this work, we investigate the dependence of the frac- 
tal dimension of sea-coast interfaces on the long-range 
correlations of synthetic landscapes. In order to study 
the invasion of the sea through the coast we consider a 
simple lattice invasion model. Each coast site is charac- 
terized by a resistance to erosion which is a function of its 
local lithology parameter and coastal configuration. Spa- 
tial long-range correlated surfaces are generated with the 
Fourier filtering method [?71I31) . which allows to control 
the nature and the strength of correlations. 

The manuscript is organized as follows. In Section [ll] 
we introduce the model and its relevant definitions. Re- 
sults for correlated and uncorrelated distributions of the 



lithology are discussed in Section III Finally, conclusions 
are draw in Section llVI 
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II. THE MODEL 

In a recent work, Sapoval et al. 21j proposed a model 
that explicitly split the dynamics of erosion into two dif- 
ferent mechanism based on their characteristic time scale. 



FIG. 1. (color online) Snapshots of typical configurations obtained with the erosion model for uncorrelated lithology on a lattice 
with 512^ sites. The sea sites are in blue (dark grey) and all the other sites are land sites (light grey). Three different regimes 
are obtained: a) subcritical, weak sea force, regime (/ < /c): the coastline is rough but not fractal, b) critical regime (/ = fc): 
the coastline is a self-similar fractal, c) supercritical, strong sea force, regime (/ > /c): the coastline is no longer self-similar, 
but rather self-afRne. 



namely, slow and rapid dynamics. While the former is 
mainly related to chemical processes, the latter is solely 
due to mechanical erosion. Here we only focus on the 
rapid dynamics, since the slow one takes place on a time 
scale which is beyond the scope of this study. 

For simplicity, the system is mapped onto a regular 
square lattice, where each site can be either a sea or 
a land site. For land coast site, a lithology parame- 
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FIG. 2. (color online) Dependence of the sea-trapping prob- 
ability, Ps, on the erosion force /, for landscapes generated 
with a uniform distribution of the lithology parameter £. A 
transition is observed from a trapped, subcritical, regime - 
for weak sea force - to a perpetual, supercritical, invasion 
regime - for strong sea force. The transition occurs for 
fc = 0.523 ± 0.001 - estimated from the crossing of the lines. 
Each curve corresponds to a different system size L^, with 
L = {128, 256, 512, 1024}, and results have been averaged over 
{800, 400, 200, 100} samples, respectively 



ter £i is assigned, which coarse-grains several geologi- 
cal properties, characterizing the mechanical interaction 
with the sea. We assume that an island site completely 




FIG. 3. Size dependence of the mass of the coastline - num- 
ber of sites - at the critical sea force, f = fc- A fractal 
dimension D = 1.34 ± 0.01 is obtained from the best fit of 
the data points. Square lattices, with uniform distribution 
of the lithology parameter, of size have been considered 
with L = {32, 64, 128, 256, 512, 1024}. Resuhs have been av- 
eraged over {3200, 1600, 800, 400, 200, 100} independent real- 
izations. The inset shows, for the yardstick method, the num- 
ber of sticks needed to cover the line as a function of the 
inverse of the stick size, for a system with L — 1024, av- 
eraged over 10^ samples. From the fit, a fractal dimension 
D = 1.33 ± 0.01 is obtained. From both methods the fractal 
dimension of the coastline in the critical regime is estimated 
to be L» = 1.33 ±0.01. 



surrounded by the sea is more fragile than a similar one in 
the coastline, so that the resistance to erosion depends 
on the local coastal configuration as, 



n = 



(1) 



where rii is the number of neighboring sea sites. Ac- 
cordingly, a coast site with only one sea neighbor has 
resistance = €i, while one completely clipped by water 
has resistance — During the invasion process, the 
sea penetrates in the coastline with a constant force /. 
Initially, the coastline, defined as the interface between 
sea and land, is a straight line at the bottom of the sys- 
tem and periodic boundary conditions are applied in the 
horizontal direction. At each iteration, all coast sites in 
the neighborhood of a sea site, and with a resistance to 
erosion below the sea force, become part of the sea. The 
proposed model corresponds to a collective invasion per- 
colation where the resistance is a function of the number 
of neighboring sea sites. For = li the resistance is 
solely dependent on the lithology parameter and ordi- 
nary percolation is recovered. The introduced weakening 
mechanism promotes the erosion of earth filaments in the 
coast, typically observed for ordinary percolation. 



III. RESULTS 

A. Uncorrelated lithology 

Let us start with the study of the model on an uncor- 
related lithology, where the lithology parameter is uni- 
formly distributed in the range < £ < 1. Figure [l] 
shows a typical configuration of the system for different 
sea forces / and obtained for the same distribution of the 
lithology parameter. Two different regimes are observed 
based on the sea force, namely, a weak and a strong 
regime. For weak sea force, as shown in Fig. [ija), corre- 
sponding to low values of /, only few sites are eroded and 
invaded by the sea. After some finite number of steps, 
the sea is completely trapped by the coastline - the resis- 
tance of all coastline sites is greater than the sea force - 
and the sea cannot invade further. For strong sea force. 
Fig. [TJc), the coastline perpetually grows, without be- 
ing ever trapped. Nevertheless, some islands remain that 
will never be destroyed, i.e., their coastline sites have a 
resistance above the sea force. As discussed before, in 
reality two different mechanisms occur: the mechanical 
and the chemical erosion. In this work, we solely account 
for the former one, however, over larger time scales, due 
to the latter mechanism the sea is permanently eroding 
the coastline. 

The transition between the weak and the strong sea 
force regimes occurs at a critical force f = fc- To deter- 
mine this critical force, we analyze the sea-trapping prob- 
ability defined as the probability, for a given sea force, for 
an erosion process to be limited to a finite number of iter- 
ations. In Fig.[2]we show this probability as a function of 
the sea force, for different system sizes. For the weak sea 
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FIG. 4. (color online) Finite-size scaling of the time depen- 
dence of the coastline roughness in the supercritical regime 
(/ = 0.9). Square lattices with Lx x Ly sites have been con- 
sidered, with Lx = L and Ly — 1281/. The time t has been 
rescaled by , where z is the dynamic exponent, and the 
roughness has been rescaled in units of 1/°", where a is the 
roughness exponent. The best data collapse is obtained for 
a = 0.48 ± 0.04 and z = 1.55 ± 0.11, a value that is con- 
sistent with the Kardar-Parisi-Zhang universality class |32) . 
The straight line has slope /3, corresponding to the growth 
exponent. The inset shows the time dependence of the coast- 
line roughness. Systems with L — {32, 64, 128, 256} have been 
sampled over {3200, 1600, 800, 400} independent runs, respec- 
tively. 



force, f < fc, since the sea invasion is always halted after 
some steps, the probability is one. On the other hand, 
in the strong sea force regime, f > fc-, the sea erosion is 
never stopped. From the crossing of all lines (different 
system sizes), it is possible to estimate the critical force 
as fc = 0.523 ± 0.001. At the critical force, the front 
of the coastline is self-similar with a fractal dimension 
D = 1.33 ± 0.01. To obtain the fractal dimension, we 
considered two different methods: the scaling of the set 
of coastal sites with system size (main plot of Fig. [3| and 
the yardstick method (inset of Fig. |3| [33] . The critical 
force obtained under the proposed weakening mechanism 
is lower than the one corresponding to ordinary percola- 
tion {fc ~ 0.593), recovered when the resistance is solely 
dependent on the lithology parameter, in the absence of 
weakening. In such case, the obtained fractal dimension 
is compatible with 7/4, the one of the hull of the perco- 
lation cluster pS] . 

For strong sea force, / > fc, the sea constantly erodes 
the coast, being never trapped, and the resulting coast- 
line front is self-affine. Moreover, our results indicate 
that the statistical properties of the sea invasion pro- 
cess belong to the same universality class of interfaces 
that obey the Kardar-Parisi-Zhang (KPZ) equation |32j . 
To characterize the main interface (coastline front), we 
study the time evolution of the roughness, W, namely. 
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FIG. 5. Size dependence of the growth exponent, p. These 
results indicate that, in the thermodynamic limit, /3 converges 
to (3 = 1/3 (dashed line) as in the Kardar-Parisi-Zhang uni- 
versality class. Systems with L = {32, 64, 128, 256} have been 
sampled over {3200, 1600, 800, 400} independent runs, respec- 
tively. 



the standard deviation of the interface width, defined as. 



W = 



(2) 



where yi is the vertical position of the interface and {y) 
its mean value over all columns. Results for different 
system sizes are shown in the inset of Fig.|4j Initially, the 
roughness is an algebraic function of time, W ^ , where 
/? is the growth exponent. At a certain crossover, t = tx, 
lateral correlations in space resulting from the invasion 
process become so large, as compared to the finite size 
L, that the roughness reaches a saturation value, Wsat 
[53] . With the system size, Wgat diverges as Wsat ~ L" 
and the crossover time scales as ~ i^, where a is 
the roughness exponent and z the dynamic exponent. 
Therefore, a complete finite-size scaling can be obtained 
with the ansatz 1351 , 



W{L,t) = L" 



(3) 



The main plot of Fig.|4]corresponds to the finite-size scal- 
ing of the roughness. From the data collapse, we obtain 
the values a = 0.48±0.04 and z = 1.55±0.11, which are 
consistent with the exponents a = 1/2 and z ^ 3/2 from 
the KPZ universality class [32l[36]. The size dependence 
of the growth exponent is shown in Fig. [5j Our results 
indicate that, in the thermodynamic limit, the growth ex- 
ponent converges to the expected KPZ value of /3 = 1/3. 
For strong sea force, when the weakening by neighboring 
sea sites is neglected, the supercritical regime of ordinary 
percolation is recovered. 



B. Correlated lithology 

In this section we discuss the case where the disordered 
coastal landscapes possess spatial long-range correlations 
in the lithology parameter i. As in previous studies [251 - 
1511 , spatial long-range correlated distributions are 

obtained with fractional Brownian motion (fBm) pi I27j. 
To achieve the intended correlated distribution, Fourier 
coefficients are generated in the reciprocal space of fre- 
quencies /, according to a power-law spectral density, 
namely. 



(4) 



where d is the spatial dimension (d = 2 in this work), 
and the inverse Fourier transform is applied to obtain 
the distribution in real space. Several samples are then 
generated and the distributions truncated between —3a 
and 3a, where a is the standard deviation. The trunca- 
tion is such that values outside this range are assigned 
to be ±3(7, keeping the original sign. Finally, the distri- 
bution is rescaled in the interval [0:1]. Each distribution 
is characterized by a Hurst exponent H - related to the 
spectral exponent by w ^ 2H + d - such that, for two di- 
mensions, the correlations are negative for < 7? < 1/2 
and positive for 1/2 < 77 < 1. The former case means 
that neighbors of a strong site are, on average, weak sites, 
whereas in the latter case they are typically strong sites. 
For H = 1/2, the classical Brownian motion is recovered 
where the increments are uncorrelated but the obtained 
lithology is still correlated. The uncorrelated distribu- 
tion of lithology is solely obtained for a constant spectral 
density, with w = and H = —d/2 {H = —1 in two 
dimensions). Further details about the adopted method- 
ology can be found, for example, in Ref. [3T] . 

Figure [6] shows snapshots of the system, at the critical 
force, for three different values of the Hurst exponent, 
namely, 0, 0.2, and 0.8. In the first two the spatial cor- 
relations are negative, while the latter corresponds to a 
strong positive correlation. By increasing the Hurst ex- 
ponent, islands become less frequent and the interface 
smooths out. In the limit of very strong positive correla- 
tions, the coastline converges to a non-fractal object with 
dimension one. In Fig. [7] we show the fractal dimension 
of the coastline front as a function of the Hurst expo- 
nent. The fractal dimension decreases continuously from 
D = 1.34±0.02, for 7J = 0, to I? = 1.04±0.03, for strong 
positive correlations. 



IV. CONCLUSIONS 

In this work, a model is introduced to study the me- 
chanical erosion of the coastline by the sea. Despite shar- 
ing some features with invasion percolation, in this model 
the update occurs over the entire interface and not only 
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FIG. 6. (color online) Snapshots of the initial distribution of lithology parameter (a, b, and c) and the corresponding eroded 
system (d, e, and f) at the critical force (/ — fc), for three landscapes generated with long-range spatial correlations, for 
different Hurst exponent: a) and d) H = 0.0; b) and e) H = 0.2; and c) and f) H = 0.8. The color scheme of coast sites 
represents the value of the lithology parameter I and sea sites were assigned with I = 0. Pictures have been obtained for square 
lattices with 512^ sites. 
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FIG. 7. Critical coastline fractal dimension D as a func- 
tion of the Hurst exponent H. A continuous decrease from 
D = 1.34 at H = to D = 1.04 at H = 1.0 is ob- 
served. The dashed line stands for the fractal dimension of 
the uncorrelated case (D = 1.33). Each point is obtained 
from the size dependence of the number of coastline sites for 
system sizes L = {32,64,128,256,512,1024} averaged over 
{3200, 1600, 800, 400, 200, 100} samples. 



the coastline and the eroding process stops after some 
time, in the latter, erosion is perpetual and belongs to 
the Kardar-Parisi-Zhang universahty class. The transi- 
tion between these two regimes occurs at a critical sea 
force, characterized by a fractal coastline front. For an 
uncorrelated distribution of the coast lithological prop- 
erties, the fractal dimension of this interface is related 
with the accessible perimeter of the percolation cluster, 
whereas for coasts with long-range correlation, the fractal 
dimension changes with the Hurst exponent. In Ref. [21] 
a model with a self-stabilization mechanism is proposed 
to explain how the system is found at criticality. Here, 
we solely focus on the effect of correlations on the prop- 
erties of the critical state without self-organization. Yet, 
this study clarifies the relation between the spatial cor- 
relations in the lithology and the multitude of fractal 
dimensions observed for coastlines. As an extension of 
this work, one might consider the effect of disorder in 
the strong sea force regime. 
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